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This article addresses the problem of two- and higher dimensional pattern matching, i.e. the 
identification of instances of a template within a larger signal space, which is a form of registration. 
Unlike traditional correlation, we aim at obtaining more selective matchings by considering more 
strict comparisons of gray-level intensity. In order to achieve fast matching, a nonlinear thresholded 
version of the fast Fourier transform is applied to a gray-level decomposition of the original 2D 
image. The potential of the method is substantiated with respect to real data involving the selective 
identification of neuronal cell bodies in gray-level images. 



Pattern correlation and image registration represent 
essential operations in statistical physics and image anal- 
ysis. Given a template reference, one is asked to find its 
related instances in a usually larger signal space. Often 
called template matching, and sometimes image registra- 
tion, such a task has usually been understood and imple- 
mented in terms of the convolution or correlation oper- 
ators (e.g. [ESHi)- Thus, given a reference pattern t, 
the identification of its instances along the space s is ob- 
tained by looking for maximum peaks of R{t, s) = t o s, 
namely the correlation between images t and s. While 
fast, previous FFT based methods of registration have 
shown significant accuracy problems. Being a bilinear 
operator, the correlation is proportionally affected by the 
intensity of gray scale images t and s. As a consequence, 
the presence in s of even a single point of high inten- 
sity can produce a false matching, substantially reducing 
the discriminatory power of the whole approach. Indeed, 
even when the patterns in s have similar height, the cor- 
relation tends to produce several significant false alarms. 
Figure ^ illustrates the rather limited power of correla- 
tion for pattern recognition. The result of the correlation 
between the template in (a) and the signal in (b) is pre- 
sented in (c). It is clear that almost any pattern within 
a broad range of intensities can produce false alarms. 
Such a poor discrimination is ultimately a consequence 
of the tolerance allowed by the correlation to perturba- 
tions of the original patterns. Though such a tolerance is 
essential to cope with the many distortions found in ex- 
perimental data, it should be kept to manageable levels 
so as not to compromise the overall discriminability. 

The current article presents a non-linear spectral ap- 
proach to template matching which provides an excel- 
lent balance between discriminability and level of toler- 
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FIG. 1: The correlation between the template in (a) and the 
signal in (b), shown in (c), allows poor discrimination between 
the patterns. The spectral matching, shown in (d), provides 
excellent discrimination of the original template. 



ance to distortions. The method involves the decom- 
position of the template and signal into binary subsig- 
nals, which are then strictly compared. In order to en- 
hance speed, the comparisons are performed by using the 
Fourier transform in order to obtain the correlation be- 
tween the quantized signals. Although a similar scheme 
has been previously applied to ID biological sequence 
matching (e.g. |j, la|)i its generalizations to images and 
higher dimensional patterns are discussed in the present 
work. 

Initially, the image signals t and s, originally with G 
gray levels, have their intensity resampled to g < G gray 



levels. This is illustrated in the diagram in Figure 12 
Note that each of the resulting subsignals ti and Si are 
binary, in the sense of being composed exclusively by val- 
ues zero or one. The strict comparison between each pair 
{ti, Si), i.e. Ci = tiQ Si is performed by using the Fourier 
transforms. Let Ti and Si be the respective Fourier trans- 
forms of ti and Si. The strict matching between ti and 
Si can be conveniently obtained as the inverse Fourier 
transform of TS. Note that the peak value of r^ will cor- 
respond to the maximum number of ones in the elemen- 
twise (Kronecker) product between ti and Si considering 
all possible lags H . The final result of the strict match- 
ing between t and s is obtained by adding all partial 
matchings, i.e. M{t,s) — '^iCi. Therefore, the choice 
of the number of quantization levels g provides the means 
for controlling the discriminative power of the method, 
with enhanced discrimination being obtained for larger 
numbers of quantization levels g. Figure ^d) shows the 
strict matching of the patterns in (a) and (b) assuming 
16 levels of binary quantization. The definite superior- 
ity of the proposed methodology can be readily verified 
by comparison with the traditional correlation n result 
shown in (c). 
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FIG. 2: Diagram presenting the main blocks in the strict pat- 
tern correlation method: the template and signal are quan- 
tized into g binary subsignals which are subsequently com- 
pared by using the Fourier transform . 



each thresholded version Si of the original image S (with 
maximum gray level G) as S; = (1 -|- sign{S > iAG))/2, 
where AG = G/{g — I); (ii) obtain the partial strict 
matchings, i.e. Ci = ti Q Si = F~^{F{ti)F{ti)), where 
F{-) is the Fourier transform; and (c) obtain the global 
strict match M = J2^=i ^i- 

The potential of the proposed methodology is here il- 
lustrated with respect to real images of neuronal cell 
bodies (somata) obtained from the public Brainmaps 
dabatase Q. Figure |2Ia) shows a portion of the Macaca 
mulatta cortex (coronal slice), containing a number of 
Nissl stained neuronal cell bodies with approximate but 
varying sizes, shapes and gray level intensities. The neu- 
ronal cell body circled in Figure Ola) has been chosen 
as the prototype (template) for the matching with the 
other cells. Note that this specific template cell presents 
some oblique elongation. The result obtained by apply- 
ing traditional correlation followed by threshold (at half 
of the maximum gray level) and peak detection is shown 
in Figure effig:exp(b). It is clear from this figure that the 
known shortcomings of the traditional correlation under- 
mine the discrimination power of the method, resulting 
in the detection of almost every cell body in a way which 
is largely irrespective to their shape or gray-level. In 
addition, some cells have been missed while other struc- 
tures have implied in overlapped peaks. The result of the 
strict matching considering the same template is shown 
in Figure ini^c). Note that only the cell bodies more sim- 
ilar to the template have in (a) in shape and gray level 
intensity have been mapped into well defined peaks in 
(b). Although remarkably good, this result can be fur- 
ther enhanced by thresholding the image in (b) at level T, 
detecting the obtained connected groups of points, and 
calculating their respective centers of mass. Figures |2Ic) 
shows the obtained centers of mass superimposed onto 
the original image (T at half of the maximum gray level 
value in (b) and g = 4). Only those cell bodies more 
similar to the template have been identified, without su- 
perposition. It is also clear from this result that, in addi- 
tion to providing a good discrimination power, the strict 
matching also caters for reasonable level of tolerance to 
pattern alterations. 

All in all, we have shown that high quality matching 
between patterns can be achieved by using a non-linear 
spectral method involving the decomposition of the sig- 
nal into subsignals followed by their strict comparison. 
While such a methodology had been previously applied 
with encouraging success for the matching between bio- 
logical sequences in bioinformatics, we have shown that 
its superiority is also manifest for the case of two- and 
higher dimensional patterns. 



Next we turn our attention to the application of the 
strict matching methodology to patterns in the plane 
{B?). The methodology, which is a direct extension of 
the ID matching strategy described above, therefore in- 
volves the following steps, given the specific number of 
quantizations levels g: (i) for i = 1,2, ...,g determine 
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FIG. 3; Image of Nissl stained neuronal cells (Macaca mullata) (a). The template is shown inside the circle. The correlation 
between this template and the whole image is shown (b), and the strict matching results are shown in (c). The respectively 
detected cells superimposed onto the original image are presented in (d). 



